Study of the K\ — > t^tx^vv decay 
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The rare decay K% — > -k°-k°vv was studied with the E391a detector at the KEK 12-GeV proton 
synchrotron. Based on 9.4 x 10 9 K% decays, an upper limit of 8.1 x 10 -7 was obtained for the 
branching fraction at 90% confidence level. We also set a limit on the K\ — >• 7r°7r°X (X — > invisible 
particles) process; the limit on the branching fraction varied from 7.0xl0~ 7 to 4.0xl0 -5 for the 
mass of X ranging from 50 MeV/c 2 to 200 MeV/c 2 . 



I. INTRODUCTION 

In the standard electroweak theory, flavor-changing 
neutral current (FCNC) processes are strongly sup- 
pressed and can only occur via higher-order effects. 
Hence, these processes will be sensitive to unknown par- 
ticles or interactions that contribute in higher-order loop 
diagrams. Such processes are ideal places to look for sig- 
nals of new physics beyond the Standard Model (SM). 
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In the SM, the FCNC K° L ir°ir°vi> decay is pre- 
dominantly a CP conserving process. Its branching frac- 
tion is sensitive to the real part of the s — > dvv tran- 
sition amplitude, while the related decays K\ — > ir^vv 
and K + — > -k^vv sense the imaginary part and absolute 
value, respectively. Like these decays, K\ — > tt tt vv 
is theoretically clean and uncertainties in the hadronic 
matrix element can be removed by using the measured 
branching fraction of its relevant semileptonic decay 
K+ Tr°ir°e + v. 

The SM predicts the branching fraction to be (1.4 ± 
0.4) x 10" 13 0, 0]. Although the prediction is solid in 
the SM, there is a possibility of enhancements from new 
physics contributions. In fact, phenomenological analy- 
ses give constraints on possible enhancements by up to 
an order of magnitude within the allowed range of new 
physics parameters from known kaon decays, including 
the measured branching fraction of K + — > tt + vv 

In addition, a new particle X which decays into in- 
visible particles, can also contribute to the three body 
decay if 9 — >• 7r°7r°A. There is also a possibility of hav- 
ing new pseudoscalar particles, as predicted by several 



FIG. 1. Schematic cross sectional view of the E391a detector. The entrance of the detector is at "0 m." 



supersymmctry models. Searches for the X — > 77 or 
X — > /i/I modes and comparisons with the models have 
been made [H, and similar test should also be tried in 
the vv final state. 

The E391a experiment at the KEK 12-GcV proton syn- 
chrotron made the first search for the -> ifiifivv de- 
cay, based on the data sample taken in the first stage 
(Run-I). An upper limit on the branching fraction of 
4.7xl(T 5 at the 90% confidence level (C.L.) for K° L -» 
-k^-k^vv was obtained Unfortunately the Run-I data 
were compromised by large neutron backgrounds from a 
vacuum membrane that hung in the beam. This article 
is based on the E391a data taken in the periods from 
February to April in 2005 (Run-II) and from November 
to December in 2005 (Run-III) , after correcting the prob- 
lem. These data had six times more K\ decays and were 
analyzed with improved methods. 



II. E391A DETECTOR 

The E391a detector wasprimarily designed to search 
for the K\ — > tt^vv decay It had a Csl calorimeter 

to detect two photons from a ir° decay and hermetic veto 
counters to ensure that no other visible particles existed 
in the final state (Fig. [T|). To avoid backgrounds from 
interactions of the beam particles with air, most of the 
detector components were placed in a vacuum vessel. 

The Csl calorimeter consisted of 496 blocks of 7 x 7 x 30 
cm 3 undoped Csl crystal. A 12 x 12 cm 2 beam hole was 
located at the center of the calorimeter to allow the beam 
particles to pass through. The main barrel (MB) and 
front barrel (FB) counters were sandwiches of lead plates 
and plastic-scintillation counters with 13.5 Xq and 17.5 
Xo, respectively, and formed cylindrical walls surround- 
ing the Kl decay volume. 

Collar shaped counters (CCOO, CC02— 07) were placed 
around the beam line for vetoing photons along the beam 
axis. Charged particles that would hit the Csl calorime- 
ter were identified and rejected by energy deposits in 
a charge veto (CV) scintillation-counter hodoscope, lo- 
cated 50 cm upstream; it covered the front face as well 



as the outer wall of the beam-hole area. 

Beam-hole charge veto (BHCV) and back-anti (BA) 
counters were located in the downstream region along 
the beam axis and outside of the vacuum vessel. The 
BHCV consisted of eight 3-mm-thick plastic-scintillation 
counters, which were arranged to fully cover the down- 
stream area of the beam hole. The BA was made of six 
supcrlayers, each consisting of a lead/scintillator sand- 
wich and quartz blocks for Run-II. For Run-III the BA 
had five superlayers, where the lead/scintillator sandwich 
part was replaced by the Lead Tungsten Oxide (PWO) 
crystals; this change was the only difference in the E391a 
detector system between the Run-II and Run-III periods. 

Data acquisition was made with a hardware trigger 
that required two or more electromagnetic showers in the 
Csl calorimeter with total energy > 60 MeV, and no ob- 
vious activity in the CV and other veto counters. Loose 
on-line veto requirements were imposed, which ensured 
flexibility of setting energy thresholds consistently for the 
real data and the detector simulation during the off-line 
analysis. The K\ — > tt tt i/D decay should produce ex- 
actly four showers in the Csl calorimeter and nothing else 
in all the other counters. These hits should also satisfy 
an in-time requirement to suppress backgrounds. 



III. DATA REDUCTION 

The experimental signature of ifiifivv is four 

photons in the final states with a non-zero missing mass 
and a transverse momentum. Photons were registered in 
the Csl calorimeter as energy deposits in adjacent blocks 
(clusters). The most probable hit position of the photon 
was obtained from the shower energy distribution among 
the Csl blocks. Events with four photon clusters in the 
Csl were selected for further off-line processing. 

Reconstruction of a pair of photons was made by as- 
suming that they were from a 7r° that decayed on the 
beam axis (z-axis). Multiple pairings of four photons to 
reconstruct two 7r°s were eliminated by taking the combi- 
nation that gave the minimum x 2 value in the difference 
in z vertices. A cut on the minimum \ 2 for the common 
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z vertex V z was effective for good tt° identification. Also, 
requiring a large difference in the lowest two \ 2 values 
was effective in reducing the combinatorial error. We re- 
quired that the common 7r°-decay vertex point V z should 
be in the decay volume 300 cm < V z < 500 cm. 

The veto conditions on the MB, CV, and collar coun- 
ters were carefully tuned. Among these counters, the MB 
played a major role in vetoing photons undetected by the 
Csl calorimeter. If too low an energy-threshold was im- 
posed for all of the MB counters, back-splash from the 
CsI/CV surface and electromagnetic shower leakage from 
the outer part of the Csl calorimeter to the downstream 
area of the MB would cause acceptance loss. Thus, we 
applied a tighter threshold to the upstream region of the 
MB to detect photons with a high efficiency, and a looser 
threshold to the downstream region of the MB to keep 
a large acceptance for the signal Q. Suppression of the 
fusion cluster (misidentification of two close-by clusters 
as one cluster) was made by using the neural network 
(NN) technique trained by the stand-alone photons and 
the fusion clusters from the —> 7r 7r°7r° Monte Carlo 
(MC) simulation sample. Further suppression was made 
on the shower shape of the cluster, and by examining the 
energy distribution and its RMS fraction Q. 
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FIG. 2. P T vs. effective mass (M^o) plots of the 7rV° 
system from Run-III MC simulations for (a) K\ — > 
(b) K° L ->■ 7T 7r°, (c) K° L -> rrW, and (d) Run-III real data. 
The statistics of the MC events in (b) and (c) are normalized 
to those of the real data. The signal region is defined to be 
inside of the triangular area in the plots. A curved line in (a) 
indicates the kinematical limit for the K decay. 



TABLE I. Reduction factors and resultant acceptances in each 
data reduction step (signal Monte Carlo). 



Cuts 



Reduction Resultant 
factor acceptance 



K" L decay in the fiducial volume 




1.00 






+ hardware trigger & four-cluster- 










candidate in the Csl calorimeter 


0.0499 


4.99 


X 


10 


+ tighter Csl fiducial 


0.616 


3.08 


X 


10 


+ rejection of fused cluster 


0.361 


1.11 


X 


10 


+ req. good cluster 


0.295 


3.28 


X 


10 


+ x" cut for two 7T° vertex matching 


0.606 


1.99 


X 


10 


+ other misc. requirements on Csl 


0.829 


1.65 


X 


10 


+ MB and upstream Vetoes 


0.621 


1.02 


X 


10 


+ downstream Vetoes 


0.853 


8.74 


X 


10 


+ req. signal box (after all the cuts) 


0.346 


3.02 


X 


10 



IV. MC SIMULATION 



A GEANT3-based [8| MC simulation was performed 
to define the signal region for K\ — > tt°tt°i/D and to es- 



timate the acceptance. The — > tt tt°i'D decay was 
generated with the V-A interaction, and the events were 
processed with the same reconstruction code as the real 
data. Among the incoming K\ particles to the detector, 
2.66% decay within the fiducial decay volume. A total 
of 1 x 10 8 incoming K\ events were generated for both 
Run-II and Run-III. The generated events were overlaid 
with accidental hits on the counters accumulated by spe- 
cial triggers during the data taking. After applying the 
selection criteria described above, a two-dimensional dis- 
tribution of the ir°ir° effective mass M^o^o and the trans- 
verse momentum Pt for the accepted events was exam- 
ined. The result of the simulation for K^ L — > ifiifivv in 
the Run-III condition is shown in Fig. [2ja). 

The signal region was defined as a right triangular 
shown in Fig. [5J The vertical edge of the triangle is 
located at the lower bound of the 7r°7r° effective mass 
(0.268 GeV/c 2 ), and spanned the Pt values between 0.1 
to 0.2 GeV/c. The horizontal edge defined the lower 
bound of Pt (0.1 GeV/c) to avoid contamination from 
the K\ — > 7r°7r°7r decay; it spanned the M^o n o values 
between 0.268 to 0.450 GeV/c 2 . 

The detector acceptance was calculated by taking the 
ratio of the number of events in the signal region to the 
number of particles that decayed in the fiducial vertex 
region. A list of the reduction factors that contributed 
to the acceptance after the K®s decayed in the fiducial 
volume are summarized in Table [U The acceptances ob- 
tained for Run-II and Run-III were consistent, and the 
result was (3.02 ± 0.11 stat .) x 10" 4 for each run period. 
To study the dominant backgrounds, MC events of the 
decay with the statistics equivalent to 11.2 
and 11.9 times the Run-II and Run-III data, and the 
Kl —> 7r°7r 7r° decay equivalent to 0.80 and 2.52 times 
the Run-II and Run-Iff data, respectively, were gener- 
ated with the same detector conditions as the real data. 
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Within these statistics, no MC events of these processes 
were found in the K\ —> ifiifivv signal region for both 
the Run-II and Run-III conditions. 

The MC results of the K° L -> 7r°7r° and K° L -> 7r 7r°7r° 
decays for the Run-III condition, with the same statis- 
tics as the real data, are shown in Figs. [2(b) and (c), 
respectively, together with the Run-III real data (Fig. 
[2(d)). Thc K l 7r ° 7r ° P lot of p T vs. M„o^o showed clear 
clustering well removed from the signal region of — >• 
tt°tt°vD (Fig. [2(b)). Events from thc K° L -> 7r°7r 7r° decay 
were spread over the vs. M^o^o plot without having 
clear boundaries (Fig. [1(c)). The data plot (Fig. [2(d)) 
outside of the masked signal region was characterized by 
the sum of Fig. [2(b) and (c), except for a small number of 
events that extended to higher Pt from the K\ — > 7r°7r° 
cluster region. These latter events are considered to be 
due to Kg regeneration in the detector, and are not sim- 
ulated in the MC. They will be discussed in Sec. VI. The 
MC study showed that most events stayed close to the 
signal region except that those for M^o^o w Mr-o were 
from K° L -> 7T 7r° and K\ — > 7r 7r°7r° decays. Thc main 
background source was iC 9 — > 7r 7r 7r . 



SENSITIVITY 



The number of if 9 particles that decayed in the fidu- 
cial vertex region was estimated by detecting if 9 — > 7r°7r° 
with the same tt° selection criteria in the Csl calorimeter 
as for Kl — >• ir°ir°vD. Comparison with the MC simula- 
tion of K\ — > 7r°7r° , with the same tt° selection criteria as 
the real data, along with overlaid accidental events, gave 
the number of K\ decays as (5.48 ± 0.09 s t a t. ±0.31 SJ/S .) x 
10 9 for Run-II and (3.88 ± 0.08 stat . ± 0.21 sys .) x 10 9 for 
Run-III, respectively 

To estimate the systematic uncertainties, the fractional 
difference between the data and the MC simulation in 
each selection criteria was examined, and the quadratic 
sum weighted by the effectiveness of each of the accep- 
tance determinations was taken. The calculated val- 
ues were 5.6% and 5.5% for Run-II and Run-III, respec- 
tively. The same values of the systematic uncertainty 
were adopted to the K\ — > tt -k°vv decay. 

Based on the acceptance for the K\ — > ifitfivv de- 
cay obtained from the MC simulation, and the number 
of Kl decays in the data-taking runs, the single event 
sensitivity (S.E.S) was defined as 



S.E.S. = 



Acceptance x Number of if 9 decays 



The S.E.S. was derived as (6.04 ± 0.24 stat . ± 0A8 sys .) x 
1CT 7 for Run-II and (8.53 ± 0.35 stat . ± 0.66 sys .) x 1CT 7 
for Run-III. The combined S.E.S. was (3.54 ± 0.1CW. ± 
0.28 sys .) x 10~ 7 . The sensitivity reported here is consid- 
erably improved from that of our previous analysis 
In addition to the six-times larger statistics, the accep- 
tance factor was improved for several selection criteria. 
Examples include the adoption of a looser cut for \ 2 m 
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FIG. 3. Kinematic boundaries of Pt vs. M^o^o for the bifur- 
cation analysis. 



TABLE II. Numbers of the events in the real data (Run-II and 
Run-III combined) selected by the combinations of two cut 
sets, A and B, at various regions defined in Fig. [3] They were 
used for the estimations of Nbkg in the bifurcation method. 
Nab is the observed number of events in those regions. 



Region 


N AB 


Nab 


Nab 


Est. N bk g 


Nab 


Signal 


3 


381 


2508 


0.46±0.26 




Low Pt 


136 


6330 


41418 


20.8±1.8 


21 


2nd low Pt 


1151 


17525 


105455 


191.3±5.9 


229 


High P T 





110 


604 


0.0 


1 


2nd high P T 


1 


5 


41 


0.12±0.13 





High Mass 


6 


283 


1086 


1.56±0.65 


12 



the two-7r° vertex matching, and optimizing the criteria 
used to select good photon clusters in the Csl calorime- 
ter, both of which were possible in a cleaner environment 
with less neutron background in Run-II and Run-III. 



VI. BACKGROUND ESTIMATION 

A "blind analysis" technique was used to minimize the 
contribution from backgrounds without examining the 
candidate events in the signal region. Among the back- 
ground sources, K\\ — \ 7r°7r° was fully reconstructed and 
rejected primarily by kinematic variables (Pt-, M^o^o) of 
the 7r 7r system. However, K\ — > 7r°7r 7r° would become 
a background if two photons were missed due to finite 
photon- veto inefficiencies in the detector. The kinematic 
cuts could not be applicable to this dominant background 
source, and suppression depended strongly on the photon 
veto conditions. 

To avoid the difficulty of having enough statistics in 
the MC simulation of the K a L — > 7r°7r 7r° background in 
the signal region, we adopted the "bifurcation method" 
to estimate the number of background events in the sig- 
nal region from the real data [![[§]• We imposed several 
selection criteria (cuts) to identify the K 9 — > Tr°Tr°i/D sig- 
nal from backgrounds, and thc roles of these cuts were 
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Proper time(s) 

FIG. 4. Proper-time distribution of the 7r°7r° system for the 
events with P T > 0.1 GeV/c and 0.45 < M n o n o < 0.55 
GeV/c 2 with a loose requirement for Vz for the Run-II and 
Run-III, assuming the production point was on the down- 
stream edge of CC02. The solid line is an expectation from 
the K% -» 7r°7r° decay. 



categorized in two groups, namely cut A and cut B. The 
cut A consisted of "veto" requirements to ensure no par- 
ticles other than from the CsTcalorimeter region, and the 
cut B consisted of the cuts to select well-reconstructed 
7r°s from pairs of photon clusters in the calorimeter. 

If the two categories of cuts are uncorrelated, the num- 
ber of events that passed both of the cuts A and B would 
be described as 

N AB = NqP(AB) £* N P(A)P(B), 

where Nq is the number of events after the basic "setup 
cuts" are imposed prior to the cuts A and B, and P(A), 
P(B) and P(AB) are the probability of passing the con- 
ditions A, B, and AB, respectively. The number of esti- 
mated background events N^kg is derived as 

N bkg = N P(A)P(B) 

_ N P(A)P{B)N P{A)P{B) 



N P(A)P(B) 
N P(AB)N P(AB) 
N P(AB) 



n ab n Ab 



N 



AB 



where the symbols A and B are the inverse logic of A 
and B, respectively. Thus N^g in the signal region after 
imposing all of the cuts A and B could be estimated from 
the number of events obtained by the combination of AB, 
AB, and AB. 

The kinematic region defined by Pr and M^o^o of the 
7r 7r system was further divided into six as shown in Fig. 
13 For each region, Nf,k g was estimated by the formula 
above. 

Table HI] summarizes the background estimates for the 
Run-II and Run-III data combined. Various kinematic 
regions were tabulated to cross-check the validity of this 
method. The numbers of events of the real data Nab 
are consistent with the estimated values N^g, and also 
agree reasonably with the MC results of dominant back- 
grounds (sum of Kl — > 7r°7r° and 7r 7r°7r° decays), after 



TABLE III. Comparison of MC events of background pro- 
cesses (Kl — > tv°tv° and Kl — > 7r c V°7r )) after imposing cuts 
A and B and real data in the kinematic regions given in ta- 
ble |TTJ The MC results are normalized to the Kl flux for the 
real data. 



Kl- 


■+ ttV Kl - 


+ ttW 


Sum of Real 


Region 


(MC) 


(MC) the two MC data 


Signal 


0.0 


0.0 


0.0 


Low Pt 


0.5 


15.6 


16.1 21 


2nd low Pt 


11.9 


203.1 


215.0 229 


High P T 


0.1 


0.0 


0.1 1 


2nd high P T 


0.0 


0.0 


0.0 


High Mass 


0.2 


0.0 


0.2 12 



imposing the cuts A and B (Table IIII[) in most of the 
regions around the signal region. 

The real data (12 events) in the "High mass" region 
exceeded the estimated number of events (1.56 ± 0.65 
events). These events are seen in Fig. Efd) with M T o w o 
around Mr-o, extending to the higher Pt region (>0.1 
GeV/c). Fig.2]shows the proper-time distribution of the 
7r°7r° system with P T > 0.1 GeV/c, 0.45 < M w o n o < 
0.55 GeV/c 2 , and a looser V z requirement for the Run- 
II and Run-III combined, assuming that the production 
point was the downstream edge of the CC02 counter. The 
proper-time distribution (Fig. |4|) is consistent with that 
of expected from the Kg —> 7r°7r° decay, implying that 
the events are from regeneration of Kg in the upstream 
materials and subsequent 7r°7r° decay. 

Because the final state of the K s — > ir°ir° process was 
fully reconstructed, these events were insensitive to the 
photon-veto cuts (cuts A) and also to the Csl cluster 
shape cuts (cuts B). Hence the bifurcation method was 
not applicable for estimating the backgrounds from this 
process. Nevertheless, these events might occur in the 
signal region by a mis-pairing of the photons from the 
7r°s. Assuming that the mis-pairing rate for K s — > it it 
decays was equal to that for Kl — > ir°ir° decays shown in 
Fig. [5{b), the contribution of events in the "High mass" 
region to the signal region was negligibly small (< 0.03 
events). 

The number of background events in the signal region 
was estimated to be 0.46±0.26 for the Run-II and Run-III 
combined (Table HE])- 

In the discussions above, we assumed that the cut sets 
A and B were uncorrelated. A measure of the correlation 
was given by the parameter e = P(A\B) — P(A\B), as 
the difference of the probability of satisfying the cut A 
for the events that passed the cut B, and for the events 
that passed the cut B. Corrections to the N^tg due to 
the correlation are given by 



C e = ex N» 



AB 



1 



Ni 



bkg 



N 



AB 



to first order of e [3,[l(3]. In the "Low Pt" region, which 
was closest to the signal region, e was (0.34±7.73) x 10~ 4 . 
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FIG. 5. Pt vs. M n o^o distribution for Run-II and Run-III 
combined. 



TABLE IV. Summary of the acceptance, S.E.S., and 90% C.L. 
upper limit on the branching fraction with several assump- 
tions of the mass of X in the K£ — > -k°-k°X (X — > invisible 
particles) decay. 

Mx accep- Br upper limit 

(MeV/c 2 ) tance S.E.S. with 90% C.L. 



50 


3.52x10' 


-4 


3.04x10" 


7 


7.0x10" 


7 


75 


3.50x10" 


-4 


3.05x10" 


7 


7.0x10" 


7 


100 


3.45x10" 


-4 


3.10x10" 


7 


7.1x10" 


7 


125 


2.91x10" 


-4 


3.67x10" 


7 


8.5x10" 


7 


150 


2.07x10" 


-4 


5.16x10" 


7 


1.2x10" 


6 


175 


9.18x10" 


-5 


1.16x10" 


6 


2.7x10" 


6 


200 


6.21x10" 


-6 


1.72x10" 


5 


4.0x10" 


5 



bined results of the Run-II and Run-III; no events are 
observed in the signal region. By using Poisson statistics 
with zero background events in the signal region, we set 
an upper limit on the branching fraction of K\ — > ir°ir°i/D 
to be 8.1 xlO -7 at the 90% C.L. The systematic uncer- 
tainty in the S.E.S. is not taken into account in the limit. 

From the same data sample, we also derived the up- 
per limit on the — > ir°ir°X (X — > invisible particles) 
decay assuming three-body phase space. Because the ac- 
ceptance depended on the mass of X, the upper limit 
was obtained as a function of the assumed mass of X as 
summarized in Table ITVl 

There are prospects of further improving the limit. A 
new experiment E14 is now in preparation at J-PARC 
[Tlj . aiming at studying the K\ — s> tt°vv decay with three 
orders-of-magnitude higher sensitivity than E391a by uti- 
lizing higher beam flux, longer beam time, and an up- 
graded detector with Csl calorimeter of higher efficiency 
and granularity. Further increases of the beam flux are 
planned in the upgrade. In the case of 2f£ — > ifiifivv, im- 
proving the photon cluster detection in the Csl calorime- 
ter is essential to improving the sensitivity. Hence we are 
hopeful to have significant progress on this process in the 
near future. 



Applying this value to the signal region, the correction 
to be made to the N\,k g was 0.01 ± 0.29 events. The 
magnitude of the correction is small compared with the 
quoted error for the N^ g , and this effect was not taken 
into account in deriving the final result. 

VII. RESULTS 

After all of these studies, the candidate events in the 
signal region were examined. Figure [5] shows the com- 
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